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EXECUTIVE SUMMARY 


A computational bulk- flow analysis for prediction of the force 
coefficients of hybrid fluid film bearings with angled orifice injection is 
presented. Past measurements on water-lubricated hybrid bearings with 
angled orifice injection have demostrated improved rotordynamic performance 
with virtual elimination of cross-coupled stiffness coefficients and nul or 
negative whirl frequency ratios. A simple analysis reveals that the fluid 
momentum exchange at the orifice discharge produces a pressure rise in the 
recess which retards the shear flow induced by journal rotation, and 
consequently, reduces cross-coupling forces. The predictions from the model 
correlate well with experimental measurements from a radial and 45° angled 
orifice injection, 5 recess water hybrid bearings (C=125 Jim) operating at 
10.2, 17.4 and 24.6 krpm and with nominal supply pressures equal 4, 5.5 and 
7 MPa. An application example for a liquid oxygen 6 recess/pad hybrid 
journal bearing shows the advantages of tangential orifice injection on the 
rotordynamic force coefficients and stability indicator for forward whirl 
motions and without performance degradation on direct stiffness and damping 
coef ficients . The computer program generated, hydrojet, extends and 
complements prior codes developed on phases I and II of the project. 
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ABSTRACT 

A computational bulk-flow analysis for prediction of ^Performance and 
force coefficients of hybrid (combined hydrostatic - hydrodynamic) bearings wit 
angled orifice injection is presented. Hybrid fluid film bearings b « a “"^ or 
cryogenic turbopumps offer reliabity with maximum operating life and optimum 
controllable rotordynamic characteristics at the lowest cost. However, fixed 
geometry hybrid bearings have limited hydrodynamic stability characteristics. 
Elements on water hybrid bearings with angled orif cc inj.c ion have 
demostrated improved rotordynamic performance with virtual eliminatio 
cross-coupled stiffness coefficients and null or negative whirl frequency 
ratios. The analysis reveals that the fluid momentum exchange at the orifice 
SiJc^ge produces a pressure rise in the recess which retards the shear flow 
induced by journal rotation, and consequently, reduces cross-coupling forces, 
ihe predictions from the model are compared with experimental measurement for a 
S- angUd orifice injection, 5 recess water hybrid bear ng operating between 
10.2 krpm to 24.6 krpm and with supply pressures from 4 to 7 MPa. The 
correlations include load and flow rates versus journal eccentricity, an ^ 
rotordynamic force coefficients at the journal centered position. An a PP^atio 
example for a liquid oxygen hybrid bearing also demonstrates the advantages 
an angled orifice injection design on the rotordynamic coefficients and 
stability indicator without performance degradation. 
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INTRODUCTION 


The importance of hybrid (combination hydrostatic and hydrodynamic) journal 
bearings as reliable support elements in cryogenic turbomachinery has steadily 
grown over the past few years. Hybrid journal bearings (HJBs) enable smaller and 
lighter turbopumps through no bearing DN life limitation and no sub-critical 
rotor operation. HJBs have durability, low friction and wear, accuracy of 
positioning, and large direct stiffness and damping force coefficients. The 
growth of an ”all-fluid-f ilm-bearing" technology for advanced and less expensive 
(per launching cost) turbopumps has required the development of analytical 
models and design tools, the testing of components, and the implementation of 
the technology (Pelfrey, 1995). 

Primary power cryogenic turbomachinery operates at high speeds and produces 
large fluid pressure rises (max. 30 MPa). These typical operating conditions 
determine the flow in the supporting fluid film bearing to be fully turbulent 
with dominance of fluid inertia and thermal transport effects. San Andres 
(1990-1995) provides bulk-flow analyses and computational programs for the 
calculation of bearing performance and rotordynamic force coefficients. These 
analyses have grown steadily in complexity and include unique features to model 

with exactness the flow in cryogenic liquid bearings. 

Measurements of bearing rotordynamic force coefficients and load performance 
are routinely performed at a high-speed Hydrostatic Bearing Test Facility (HBTF) 
(Childs and Hale, 1994). Tests have been conducted with water for more than 30 
hybrid journal bearings and damper seals with rotational speeds ranging from 10 
to 25 krpm and pressure differentials from 4 to 7 MPa. The facility accomodates 
state of the art instrumentation with remotely controlled testing, and includes 
an efficient real-time parameter identification method based on frequency domain 
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techniques. Kurtin et al. (1993), Franchek et al. (1994, 95), Mosher and Childs 
(1995), and Yang et al. (1995) report extensive experimental data for the static 
performance characteristics of a 5 recess HJB for the operating conditions noted 
and three different bearing clearances (76 to 127 pm). These studies demonstrate 
that predictions from the computational bulk-flow models correlate favorably 
with the experimental results. Accurate predictions depend greatly on the 
knowledge of the bearing operating clearances, and most importantly, on the 
orifice discharge coefficients. The references cited along with San Andres 
(1995a) also discuss the sensitivity of the computed predictions to variations 
in the input empirical parameters. 

Despite the many advantageous features offered by HJBs, hydrodynamic and 
"pneumatic-hammer" stability limits and two-phase flow operation are of primary 
concern for high speed operation with large pressure differentials. Fluid 
vaporization is possible since the cryogenic liquid enters the bearing (or seal) 
at conditions close to its saturation temperature. The large mechanical energy 
dissipated by the fluid motion is convected by the fluid and changes greatly its 
thermophysical state while the pressure drops to the exit plane of the bearing. 
"Pneumatic hammer" effects are avoided by appropriate selection of the flow 
restrictor, by designing bearing recesses with small volumes, and by restricting 
bearing operation to flow conditions where the pressure differential is a small 
fraction of the liquid bulk modulus (Redecliff and Vohr, 1969). This last 
condition is difficult to achieve in bearings employing liquid hydrogen due to 
its large compressibility. 

The stability of a simple rotor-bearing system is defined by its threshold 
speed and the whirl frequency ratio (WFR) . The instability of "hydrodynamic" 
type is solely due to the effect of journal rotational speed on the bearing flow 
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field. The threshold speed corresponds to the rotor speed at which a bearing is 
deprived from its effective damping and any small perturbation from an 
equilibrium position will determine unbounded rotor motions. The VFR denotes the 
ratio between the onset whirl frequency (typically the system first critical 
speed) and the threshold speed of instability. This stability indicator is 
independent of the flexibility of the rotating shaft. Plain journal bearings 
show a VFR equal to 0.50 for small to moderate operating eccentricities (light 
loads), and thus demonstrate instability at a rotational speed equal to twice 
the system first critical speed. Measurements in hybrid bearings verify closely 
the theoretical WFR prediction. In some circumstances the WFR even increases 
above 0.50, in particular for low rotational speeds and large supply pressures 

(Franchek, 1992, Franchek et al. 1995). 

The WFR=0 . 50 condition limits severely the application of HJBs to high speed, 
light weight turbomachinery, and thus, the research has concentrated on 
conceiving hybrid bearings with improved stability and without loss in centering 
stiffness and damping ability. Some of the technological advances have been the 
natural outcome of analysis and engineering design, while others follow 
empirical evidence and past experience when a mathematical model is yet to be 
crafted. Other recommended fixes to improve the hydrodynamic stability of hybrid 
bearings by reducing or eliminating the WFR are the following: 
o Use of machine roughened bearing surfaces to decrease the cross-coupled 

stiffness coefficients. Test results show a rough knurled-pat tern HJB to have 
a WFR as low as 0.30 but with a reduced load capacity and direct stiffness 
when compared to a smooth surface HJB (Franchek, 1992). 
o Use of circumferentially asymmetric pad bearings and recesses to produce 
enough anisotropy on the rotordynamic force coefficients. Measurements and 
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analysis for an engineered two pad HJB validated the concept (San Andres, 

1995b). However, this bearing configuration is highly sensitive to the 
direction of applied static loads. 

o Use of flexure-pivot, tilting pad HJBs or compliant surface (foil) journal 
bearings due to their inherent stability. San Andres (1995c, 1994) discusses 
at length these concepts and evaluates their potential for cryogenic uses. 
Flexure-pivot HJBs constitute a novel altenative and full-scale testing is 
planned for the first semester of 1996. Foil bearings have also demonstrated 
their performance in cryogenic turbomachinery (Genge et al., 1993). The 
current foil bearing technology allows only for specific loads applicable to 
secondary power cryogenic turbopumps. The interested reader should recall the 
cited references for further details. 

o Use of hybrid bearings with angled liquid injection opposing journal rotation 
to reduce the development of the circumferential flow velocity and with 
virtual elimination of cross-coupled stiffness coefficients. This concept has 
lacked firm theoretical modeling though it has proved successful in some 
applications (Tondl, 1967, Brown and Hart, 1986). Experimental measurements 
for a 5 recess HJB demonstrate that angled injection aids in reducing the 
whirl frequency ratio without decreasing the bearing centering stiffness and 
load capacity (Franchek, 1992, Franchek and Childs, 1995). 

This report presents the thermohydrodynamic analysis of real properties, 
hybrid bearings with angled orifice injection. The objective is to advance a 
computational model able to predict reliably the performance of angled injection 
HJBs in lieu of their favorable (measured) rotordynamic performance. The motion 
of a fluid thpough the thin film lands is governed by mass, momentum and energy 
transport equations for the bulk-flow velocities, pressure and temperature, 
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along with thermophysical state equations for evaluation of the cryogen material 
properties. The turbulent bulk-flow is modeled with simple friction coefficients 
and include effective film depths to accomodate for macroscopic surface 
roughness. A simple analysis for the angled injection - orifice flow reveals 
that the fluid momentum exchange produces a pressure rise in the recess which 
retards the shear flow induced by journal rotation. Zeroth-order equations 
describe the fluid flow field for a journal static equilibrium position, while 
first-order linear equations govern the fluid flow for small amplitude journal 
center translational motions. Solution to the zeroth-order flow equations 
provides the bearing flow rate, shear torque and load capacity. Solutions to the 
first-order equations determine the linearized rotordynamic stiffness, damping 
and inertia force coefficients as functions of a whirl frequency. The numerical 
predictions from the model are correlated extensively with the experimental data 
of Franchek (1992). 


ANALYSIS 

Figure 1 shows the geometry of a hybrid (combination hydrostatic/hydrodynamic) 
journal bearing and the relevant nomenclature. A liquid at high pressure (Pg) 
and inlet temperature (T g ) is supplied (radially or angled) through orifice 
restrictors and impinges into the bearing recesses with a mean pressure ( P r ) - 
The pressure field within the recesses is determined from flow continuity with 
the film lands, momentum exchange at the orifice plane and a viscous rise due to 
journal rotation. At the recess edges, an inertial pressure drop also occurs due 
to the sudden transition from the recess of depth (H r ) into the film land 
regions of small thickness (H). Past the recesses, the liquid then flows through 
the film lands and the pressure drops to the discharge value (P a ). 
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Equations of flov on the bearing film lands 

On the thin film lands flow turbulence, fluid inertia and compressibility 
effects are important. The model then assumes a fully developed turbulent 
bulk-flow of a fluid whose material properties depend on its local 
thermophysical state of pressure and temperature. The equations of mass, axial 
and circumferential momentum, and adiabatic-flow energy transport for the 
bulk-flow velocities, pressure and temperature in the bearing film lands are 
given as (Yang et al. , 1995, Kleynhans and Childs, 1995): 
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Please refer to the Nomenclature for a description of all variables. 

K =K =(k +k )/2 are the wall shear stress parameters determined as local 
y x v J B 

functions of turbulent friction factors which depend on the bearing and journal 
surface conditions and the flow Reynolds numbers relative to the rotating (Rj) 
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and stationary (R 0 ) surfaces, i.e Kj-fyRj, K B = f B' R B (Hirs ’ 1973) * The cry0genic 
liquid properties are extracted from the Benedict-Veb-Rubin equation of state as 

given in the standard data base of McCarty (1986). 

The fluid pressure at the sides of the bearing (y=±L/2) equals the discharge 
or ambient value (P ). At the interface with the bearing recesses, continuity of 

cl 

flow and pressure must be attained as detailed below. 

Angled Injection - Recess Flow and Pressure equations 

Figure 2a depicts a hydrostatic bearing recess (or pocket) with axial length 
(1) and circumferential extent (b). The figure shows the direction of the 
journal surface speed (Q*R) , and relative to this velocity the recess is divided 
into upstream (u) and downstream (d) regions. The fluid supply orifice port with 
injection angle (5) is located at a distance b u from the upstream recess edge. 
The orifice has an effective area A q normal to the feed speed V Q . Radial fluid 
supply is indicated by 5=0 while a tangential feed opposite to journal rotation 

is given by S=it/2 (90°). 

Conventional analysis of hydrostatic bearings do not calculate the flow 
field within the recess since these are typically deep and enclose large nearly 
stagnant fluid volumes. Analysis then accounts only for flow continuity with the 
film lands and determines a (uniform) recess pressure using a simple orifice 
equation based on Bernoulli's principle. The complexity of the flow field in 
hydrostatic pockets has been discussed by Hill et al. (1995) and Braun et al. 
(1993, 1995) with the aid of two-dimensional computational fluid mechanics 
analyses. Numerical results reveal the generation of hydrodynamic pressures 
within the pocket and followed by sharp inertial pressure drops at the recess 
edges. This field of study is of utmost importance for the development of a 
mature technology on hybrid bearings for cryogenic applications. 
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The analysis of angled injection - hydrostatic pockets follows here a 
simplified approach which intends to be of practical use without resorting to 
computationally intensive three dimensional flow calculations. The flow model is 
evidently crude yet it grasps the fundamental mechanisms of pressure generation 
within the bearing pockets. The favorable correlation with hybrid bearing 
experimental performance characteritics given later justifies the method used. 

A mass conservation equation at each bearing recess of area (l'b) and depth 
H is defined by the global balance between the mass flow through the orifice 
restrictor (Q ro )> the mass flow into the film lands and the time rate of change 
of liquid mass within the recess and supply line volume (V r ), i.e. 


Q = p A V 
ro r o o 


p H U.n 


dr + — (p v ) 

it 


(5) 


for r=l , 2. . . , Nrecess 


where A =C.nd^/4 is the effective orifice area with C, as an empirical discharge 
o d o a 

coefficient. T r denotes the closure of the recess with the film lands and has a 
normal n along the boundary line. At the orifice discharge plane, the mean 
recess pressure is denoted by (see Figure 2b) and given from Bernoulli's 
equation as: 


<P s -P r ) = (1/2) p r V* (6) 

Computational fluid mechanics analysis reveals that the axial pressure within 
the recess is (to a first approximation) practically uniform. Hence, modelling 
of the flow in the pocket as a one-dimensional bulk-flow bearing determines that 
the pressure difference (downstream - upstream) on a recesses is given by two 
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contributions: 


a) a viscous pressure rise (AP rv ) due to shear flow (San Andres, 1992): 



r 

b 

Q-R 

AP 

rv 

P, - P 
d u 

= u K — 0 

r xr u 2 
v H 1 
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b) a pressure drop (AP rm ) at the orifice injection plane and due to the exchange 
of fluid momentum, and simply stated as: 
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( 8 ) 


where the orifice equation (6) has been used on the right hand side of eqn (8). 
Note that for radial injection ( 5 = 0 ) there is no momentum pressure drop at the 
supply port. The viscous pressure rise depends greatly on the journal speed and 
the mean recess circumferential flow speed. On the other hand, the momentum 
exchange pressure drop is the largest for large pressure differentials ( p s _p a ) 
and tangential injection (5=90°). For simplicity the pressure field within the 
hydrostatic pocket is then taken as linear and combines the two pressure 
differences as shown pictorially in Figure 2b. Note that this simplification 
avoids the calculation of the complex flow field on the entire bearing recess. 

Finally, the entrance pressures (P e ) to the film lands in the circumferential 
(upstream and downstream) and axial directions are given by (San Andres, 1992): 


P 

e 






(9. a) 


(9.b) 
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The analysis generalizes equations (9) for uneven empirical entrance loss 
factors £ in the upstream(u) and downstream(d) sections of the recess. Equations 
(9) are used only when fluid flows from the recess towards the film lands. 
Perturbation Analysis 

Consider the motion of the journal as the superposition of small amplitude 
periodic motions of frequency («) around a static equilibrium position. That is, 
the journal center displacements are given as 


e x (t) = 


'Xo 


+ Ae, 


i»t 


* 


'Yo 


Ae, 


loot 


i=4-l (10) 


The magnitudes of the dynamic perturbations in journal displacements are small, 

i.e., | (Ae ,Ae }| «< C. The film thickness (H) can then be regarded as the 
X Y 

superposition of a steady-state (H q ) and dynamic components given by the real 
part of the following expression: 


H = H q + { Ae x h x + Ae y h ? } e iwt ( n ) 

where H q = C(y) + e Xo h x + e Yo h y ; and h x = cos(9), h y = sin(9) 

The flow field variables (U.U ,P,T), as well as the fluid properties (p,u) 

x y 

and the shear parameters (K.K) are also formulated as the superposition of 

x y 

zeroth-order and first-order complex fields describing an equilibrium for 
steady-state flow, and the perturbed condition for small amplitude dynamic 
journal motions, respectively. In general, these fields are expressed as: 


▼ = Y o + *xV*T T Y 


¥ 


itot 


( 12 ) 


Substitution of equations (11) and (12) into the flow equations (1-9) renders 
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zeroth- and first-order equations for determination of the steady-state and 
perturbed flow-fields. These equations are not reproduced here for brevity but 
can be found in their full extent in the reference of San Andres (1993). The 
bearing static and dynamic force characteristcis are evaluated once a solution 
to the flow equations is obtained. Fluid film forces (F^Fy) force 
coefficients (stiffness damping and inertia M^) are calculated by 

integration of the pressure fields over the journal surface. The appropriate 
formulae are: 


F 

a 


X- 2n 

P h R*d9*dy; 
o a 

Vo 


a=X, Y 


(13) 


R M fl +iwC « 
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P_ h R-d9*dy; 
P a 


a, 3=X,Y 


' 0 J ( 


(14) 


Numerical method of solution 

The control-volume method of Launder and Leschziner (1978) is used to 
solve the differential equations of motion. Staggered grids containing control 
volumes for the primitive flow variables (circumferential and axial velocity, 
pressure and temperature) cover the flow domain. Algebraic difference equations 
are derived on each control volume for the conservation of mass, axial and 
circumferential momentum, and balance of energy. A pressure correction equation 
is derived using the SIMPLEC procedure of Van Doormaal and Raithby (1984). A 
Newton-Raphson scheme is also used for satisfaction of the recess mass flow 
constraint. Full descriptions on the accuracy and parameter sensitivity of the 
method as applied to hybrid bearings and annular pressure seals are given in 
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past publications (San Andres, 1990-1995). The interested reader should consult 
the cited references for a detailed exposition of the numerical method used. 

The computer program generated in fortran-77 is named hydrojet and it is 
based on the original code hydrosealt developed earlier (San Andres, 1993). 
hydrojet is fully compatible with its predecessors including hydroflext and 
keeps the same basic structure and user friendly interface, hydrojet is 
available from Texas A&M University Technology Licensing Office. 


COMPARISONS TO TEST RESULTS FROM A WATER 5-RECESS HYBRID BEARING 

Franchek (1992) presents an experimental study of five hybrid bearings with 
distinctive geometrical configurations. These are namely, smooth bearings with 
radial injection and rectangular (baseline), triangular and circular recesses, a 
knurled rough-surface bearing with rectangular recesses, and a smooth surface 
bearing with rectangular recesses and a 45° angled orifice injection. The tests 
consisted of the measurement of load vs. journal eccentricity, torque and flow 
rate, and the identification of rotordynamic force coefficients (stiffness, 
damping and inertia) using a random frequency shaker excitation system. Childs 
and Hale (1994) provide a full description of the test apparatus and the 
experimental procedure. The nominal test conditions include: 


(a) 3 rotational speeds: 10.0, 17.4 and 24.6 krpm . 

(b) 3 supply pressures: 4.0, 4.5 and 7.0 MPa (600, 800 and 1000 psig) 

(c) 6 journal eccentricity ratios (e/c): 0.0 to 0.5 
at a supply temperature of 55°C (130°F). 


Franchek and Childs (1994) and Franchek et al. (1995) briefly report the 
measurements in the archival literature with comparisons to predictions from the 
present analysis for the baseline bearing (radial injection). Table 1 describes 
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the geometry of the test hybrid bearings with radial and angled (45°) injection. 
At the journal centered position the measured data for flow rate, supply and 
average recess pressures and operating clearance is also given. From these 
values, empirical orifice loss coefficients (Cj) are estimated for each test 
condition and used in all computations including journal off-centered 
operations. The values of circumferential (Re c ) and axial flow Reynolds (Re a ) 
numbers demonstrate the character of the flow within the test bearings. The 
experimental measurements along with predictions from the current numerical 
model follow. Angled and radial injection results are presented at the top and 
bottom of the figures, respectively, and with measurements depicted in broken 
lines and predicted values in continuous lines. 

Static Performance Characteristics of test bearings 

Figure 3 depicts the journal eccentricity versus applied load (W x ) for both 
bearings at a nominal supply pressure of 7.0 MPa and three rotational speeds. 

The journal eccentricity increases almost linearly with the applied load which 
is typical of externally pressurized bearings. This also demonstrates that the 
stiffness coefficients for the bearings will not vary (greatly) with the journal 
center position. An increment in the operating speed produces a rise in the load 
capacity (smaller eccentricities) with both bearings having similar load 
capacities. The model calculated results agree best with the angled injection 
tests. Figure 4 shows the predicted journal center locii for both bearings. Test 
results are not shown since these are not included on Franchek's work (1992). 
Note that for the radial bearing, the journal locus presents a greater attitude 
angle as the journal speed increases. On the other hand, the angled injection 
results show a negative attitude angle at the lowest speed (10.2 krpm) , an 
almost null angle at the medium speed (17.4 krpm), and increasing as the speed 
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reaches its highest value (24.6 krpm). The predictions reveal the fundamental 
differences between the two bearings. Thus engineering design could lead to a 
hybrid bearing free of cross-coupling effects with the appropriate combination 
of injection angle, supply pressure and operating speed. 

Figures 5 and 6 show the predicted (dimensionless) film pressures at the 
midplane of the bearing for the centered condition (no load) and at a load (V^) 
equal to 10 kN and directed towards the middle of a recess. Note that for 
concentric operation the angled injection bearing shows a recess pressure field 
significantly different from that for the radial injection HJB. As modeled, the 
exchange of momentum in the recess produces a pressure rise in the direction 
opposite to journal rotation which retards the development of the cicrcumfer - 
ential flow velocity. The effect is most pronounced at the lowest rotational 
speed (10.2 krpm). The recess pressures increase with the rotational speed 
denoting a rise in the flow resistance within the bearing film lands. 

Figure 7 depicts the recess pressure ratios {(P r - p a )/(P g - p a )} at t ^ ie j° urna l 
concentric position for the three nominal supply pressures and speeds. Recess 
pressure ratios rise with the journal speed and decrease with supply pressures 
since land flow resistance and turbulence are greater. The correlations with the 
model predictions are regarded as satisfactory except at the largest speed and 
lowest supply pressure. Note that the comparisons have been made with averaged 
test recesses pressures. Actual measured recess pressures within the bearing 
vary as much as 13% from the calculated test average. 

Figure 8 shows the predicted flow rates to agree closely with the 
measurements for both bearings at the centered journal position. The flow rate 
increases with supply pressure and decreases significantly with speed for the 
radial (baseline) bearing. On the other hand, the angled (45°) injection bearing 
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shows a less pronounced drop as speed rises while the model predictions show the 
opposite behavior. The rationale for the discrepancy is attributed to the larger 

predicted recess pressures with the model. 

Figure 9 shows the comparisons between measured and calculated drag shear 
torques at the journal centered position. It is noted that the measurements 
refer to the torque on the bearing surface which may differ from that on the 
journal. Furthermore, the measurements at the lowest speed (10.2 krpm) are 
unusually low considering the size and operating conditions of the bearing. It 
appears as if the test torque is shifted by a constant amount for both bearings. 
Nevertheless, the measurements (and predictions) show that the drag torque is 
practically independent of the external supply pressure (P g ) and increases in a 
way not (linearly) proportional to journal speed due to the turbulent flow 
character of the test bearing. Most importantly, the angled injection bearing 
presents a smaller torque than the radial injection bearing and thus, it 
indicates that shear power losses are also smaller for this hybrid bearing 
configuration. 

Dynamic Performance Characteristics at journal centered position 

The test results and numerical predictions demonstrate that the rotordynamic 
force coefficients are practically insensitive to the applied load for journal 
eccentricities to 50% of the bearing clearance. Hence, in the following, only 
force coefficients at the concentric position are presented and thouroughly 
discussed . 

The whirl frequency ratio (WFR) , a stability parameter of paramount 
importance for the application of hybrid bearings to high speed applications, is 
depicted in Figure 10. The radial bearing presents a WFR close to 0.50 for most 
operating conditions and indicates its relatively poor stability, in particular 


16 


at low speeds and high supply pressures. On the other hand, the angled (45 ) 
injection bearing shows a (large) negative VFR at the lowest speed and raising 
to the 0.5 limit as the rotational speed increases. The numerical predictions 
agree well with the measurements at the middle and high speeds, i.e. 17.4 and 
24.6 krpm. 

Figure 11 depicts the cross-coupled stiffness coefficients as 

the journal speed increases for three nominal supply pressures. In the radial 
injection bearing, the cross-coupled stiffness are always positive and increase 
with the journal speed. On the other hand, the angled injection bearing presents 
negative cross-coupled coefficient at the lowest speed. From a rotordynamics 
point of view this is a desired occurrence since then these dynamic coefficients 
provide forces damping the development of forward whirl motions. The identified 
experimental results show some variations between and -K^ which are most 
pronounced for the radial injection bearing. The numerical predictions show the 
same trends as the measurements but do not agree well with the measurements. The 
discrepancies (even with the radial injection bearing) could be attributed to 
the limitations of the bulk-flow model on handling reverse flow conditions. 

Figure 12 presents the direct stiffness coefficients < K xx = K yY^ versus the 
journal speed and nominal supply pressures. The angled (45°) injection bearing 
has larger direct stiffnesses at the lowest speed and largest supply pressure. 
However, both bearings have similar stiffness values at the largest speed 
tested. The experimental results show significant discrepancies between and 

and attributed to minute differences in the diameters of the feeding 
orifices. The numerical predictions agree very well with the measurements except 
at the lowest speed and highest pressure where the tests show an unexpected 
behavior. 
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The direct damping coefficients (C^C^) are shown in Figure 13. The model 
predictions show similar damping values for both bearings and these coefficients 
raise with the journal speed. On the other hand, the test results show damping 
to decrease with speed for the radial bearing and the opposite effect for the 
angled injection bearing. Note also that the dicrepancies between test C n and 
are significant for the radial injection bearing. In general, the model 
calculations agree best with the angled injection HJB test results, although 
direct damping is underpredicted by as much as 25% at 24.6 krpm. 

Figure 14 shows the predicted cross-coupled damping coefficients (Cyy=-Cyj) 
to increase with journal speed and with little influence of the external supply 
pressure. On the other hand, the test results show a different behavior with 
cross-damping coefficients being the largest at the middle test speed (17.4 
krpm). No conclusive remark can be made in regard to the correlation of 
prediction and identified test coefficients. 

Figures 15 and 16 show the direct inertia ( ^ yy ~^yy ^ cross-coupled 

inertia ( H ) force coefficients, respectively. The numerical predictions 

Ia A I 

demonstrate these added mass coefficients to be practically independent of 
external supply pressure and with a slow variation as the journal speed 
increases. The experimentally identified inertia coefficients are of the same 
order of magnitude as the predictions but present unique features yet to be 
fully understood. The test direct inertia coefficients do not follow a clear 
trend and show the largest values at the middle test speed (17.4 krpm). For the 
radial (45°) injection bearing, negative direct inertia coeficients at 10.4 krpm 
actually may indicate a dynamic stiffnening of the hybrid bearing. The predict- 
ions evidently do not agree well with the test results. However, Franchek and 
Childs (1994) indicate the test inertia coefficients have average uncertainties 
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of 53%, thus explaining the erratic behavior of these parameters. 

HYBRID BEARING EXAMPLE FOR A LIQUID OXYGEN APPLICATION 

Table 2 contains design data for a liquid oxygen (T g =90°K) , 6 pad/recess 
hybrid bearing operating at 25 krpm and with a pressure drop (P g -P a ) equal to 
17.9 MPa. The application corresponds to an Advanced Launching System (ALS) 
turbopump configuration (San Andres, 1995c). The analysis for journal centered 
operation investigates the effects of the angle of fluid injection on the 
performance characteristics of the bearing. An injection angle equal to 0° 
indicates radial fluid supply while an angle equal to 90° denotes tangential 
injection against journal rotation. 

Figure 17 presents the bearing whirl frequency ratio (WFR) , flow rate, mean 

recess pressure ratio (Pj-V^W and maximum P ressure <P r max) within th ® 
recess versus increasing values of the injection angle angle (6). The flow rate 
and recess pressure ratio remain practically constant for all injection angles 
while the maximum recess pressure retarding the development of journal rotation 
steadily increases and becomes a maximum for tangential injection. The most 
important result concerns the whirl frequency ratio which decreases from a value 
close to 0.60 for radial injection to approximately -0.53 for tangential 
injection. 

Figure 18 depicts the drag shear torque on the journal surface, and the 
maximum exit and mean exit fluid (T | y-±L/2 -T s > temperature differences. The 
shear torque decreases dramatically due to the effect of the fluid injection 
opposing the journal rotation. The temperatures calculated correspond to a 
thermal flow model with adiabatic journal and bearing surfaces. 

Figures 19 and 20 show the stiffness ^ X ,Y and damping < C a0>«,p-Z,T 
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coefficients versus increasing values of the angle of injection. The subindex 
(o) on the stiffness coefficients denotes values at zero frequency («=0), while 
the others are evaluated at a synchronous frequency, i.e. they contain the 
inertia coefficients and could be thought as dynamic stiffnesses equal to 

K -»^M . Note that direct stiffness and damping coefficients vary little with 

afio a0 

the angle of fluid injection. On the other hand, the cross-coupled coefficients 
decrease steadily as the orifice supply angle increases with minimum values for 
tangential injection. For injection angles greater than 25° the WFR is zero and 
then negative indicating a bearing with unlimited stability for forward whirl 
motions. 

CONCLUSIONS 

The growth of an "all-fluid-film-bearing" technology for support of 
advanced cryogenic turbopumps demands the development of models and design 
tools, the testing of components, and the implementation of the technology on 
actual hardware. Conventional hybrid fluid film bearings have demonstrated 
adequate load support, direct stiffness and damping, but suffer from limited 
hydrodynamic stability which deters their use for high speed applications and 
flexible rotating structural systems. On the other hand, experiments on hybrid 
bearings with angled orifice injection have shown virtual elimination of 
cross-coupled stiffness coefficients and null or negative whirl frequency ratios 
(Franchek, 1992). No firm analysis was available at the time of the 
measurements, and hence, further technological developments since then were 
prevented. 

The computational bulk-flow analysis for prediction of the performance 
and force coefficients of hybrid (combined hydrostatic - hydrodynamic) bearings 
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with angled orifice injection is the subject of the present research. The motion 
of a fluid through the thin film lands is governed by mass, momentum and energy 
transport equations for the flow velocities, pressure and temperature, along 
with state equations for the cryogen material properties. A simple model for the 
angled injection - orifice flow reveals that the fluid momentum exchange at the 
orifice discharge produces a pressure rise in the recess which retards the shear 
flow induced by journal rotation, and consequently, reduces cross - coupling 
forces. Zeroth-order and first-order diferential equations describe the fluid 
flow for a journal static position and dynamic perturbations, respectively 
Solution to the zeroth-order equations provides the bearing flow rate, shear 
torque and load capacity. Solutions to the first-order equations determine the 
linearized rotordynamic force coefficients at a specified whirl frequency. 

The predictions from the model are compared with experimental measurements 
for two hybrid bearings with radial and 45° angled orifice injection. The test 
bearing nominal clearance is 125 ym and operates with water at 10.2 krpm to 24.6 
krpm and with supply pressures from 4 to 7 MPa. Comparison of experiments and 
model calculations for load, flow rate and recess pressures are good and verify 
the soundness of the bulk-flow model. Correlations of model and test direct 
stiffness and damping coefficients are also favorable. The predictions show the 
same trends as the test values for the whirl frequency ratio and cross-coupled 
stiffness coefficients but large differences are apparent. Inertia force 
coefficients do not agree with the identified experimental values perhaps due to 
the large uncertainty in the measured coefficients. The experiments as well as 
the measurements demonstrate that the advantages of angled injection in hybrid 
bearings are lost as the journal speed increases and brings dominance of 
hydrodynamic over hydrostatic effects. 
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An application example for a liquid oxygen 6-recess hybrid bearing also 
verifies that a tangential angled orifice injection produces the lowest 
(negative) whirl frequency ratio and induces the largest cross-coupled 
stiffnesses which retard the development of forward whirl journal motions. 
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Table 1. Description of water lubricated orifice compensated hybrid bearing 
tested by Franchek and Childs (1994). 

No of recesses (Nrec) 5 

Clearance nominal (C) 125.4 ym (0.005 in) 

Diameter (D) 76.2 mm (3 in) 

Length (L) 76.2 mm (3 in) 

Land roughness (peak-peak) 0.33 ym (13 yin) 

Recess dimensions: square (1) 27 mm.x (b) 27 mm x 254 ym (depth) 

Supply volume Vs = 0.1289 dm 

Orifice at midplane of recess. 

diameter d=2.49 mm, C^varies (see below) 

Radial and angled injection (45°) 

Lubricant: water at Ts=328.3 °K 
Viscosity (y) 0.4929E-3 P^»s 

Density ( p) 986.26 kg/m 

Discharge pressure Pa: 0.0 MPa ( 0 psig) 


Empirical parameters: 

Entrance loss factors 1^,1 0 . 0 , Inlet swirl <x=0.5 


Test Conditions and estimated parameters at centered operation 


Angled injection: 
Speed P C 
KcDtn MPa ym 

Q 

lt/min 

Pr 

M K* 

C d 

Re 

(pflftc/y) 

Re 

(Q?nDy) 

10.2 

4.133 

5.519 

6.877 

122.8 

124.4 

124.9 

79.94 

92.21 

102.29 

1.128 

1.252 

1.434 

0.700 

0.680 

0.660 

9,998.0 

10,129.1 

10.169.8 

11,136.0 

12,845.6 

14,249.8 

17.4 

4.154 

5.521 

6.846 

120.5 

121.5 
122.7 

79.60 

91.67 

101.75 

1.866 

2.148 

2.316 

0.800 

0.763 

0.726 

16,737.3 

16,876.2 

17,042.0 

11 , 088 . 3 
11,088.9 
14,174.0 

24.6 

4.135 

5.532 

6.844 

119.4 

120.8 

117.1 

78.07 

92.21 

101.38 

2.424 

2.870 

3.206 

0.907 

0.859 

0.808 

23,447.2 

23,722.0 

23,000.0 

10 . 87 5.8 
12,942.4 
14,123.0 

Radial injection: 
Speed P C 
KcDm MPa ym 

Q 

lt/min 

PR 

Mir 

C d 

Re 

(ps£c/y) 

Re 

( p8/ JiDy) 

10.2 

4.120 

5.519 

6.889 

119.6 

119.4 

118.2 

84.89 

98.70 

110.86 

1.759 

2.163 

2.342 

0.840 

0.820 

0.790 

9.738.2 
9,722.0 

9.624.3 

11.825.0 

13.749.0 

15.443.0 

17.4 

4.118 

5.494 

6.889 

116.5 

116.5 

115.7 

82.37 

97.91 

109.84 

2.199 

2.675 

3.100 

0.904 

0.887 

0.856 

16,181.7 

16,181.7 

16,070.6 

11,474.1 

13,639.0 

15,300.7 

24.6 

4.141 

5.494 

6.904 

111.9 

115.8 

114.2 

73.11 

91.64 

105.88 

2.550 

3.164 

3.693 

0.881 

0.912 

0.900 

21,974.3 

22,740.2 

22,426.0 

10, 184.2 
12,765.5 
14,749.1 


(*) Cd values estimated from measured flow rate and average recess pressures 


Table 2. Description of 6 pad/recess liquid oxygen hybrid bearing with 
angled injection 

No of recesses (Nrec) 6 

Clearance nominal (C) 76.2 ym (0.003 in) 

Diameter (D) 92.7 mm (3.65 in) 

Pad arc length 40.45 mm (1.59 in) [50 1 

Length (L) 37.1 mm (1.46 in) 

Land roughness (peak-peak) smooth 

Recess dimensions: square (1) 192mm x (b) 19 mm x 228 ym (depth) 
Supply volume Vs = 0.0 dm 


Orifice at midplane of recess, 
diameter d=2.328 mm, C.=1.0 
Angle of injection (-10° to 90°) 


Rotational speed: 25,000 rpm (2,618 rad/s) 

Lubricant: liquid oxygen at Ts= 90 °K (supercritical conditions) 
Supply pressure P =26.71 MPa (3,874 psi) 

Exit pressure P^= 8.81 MPa (1,278 psi) 

SuddIv Exit 


0.2459E-3 

1,192 


0.2125E-3 Pa. s, 
1,160 kg/m" 


Empirical parameters: 

Entrance loss factors S^^O.O, Inlet swirl o=0.5, C d =1.0 
Circumferential flow Reynolds number, Re c =p2Rc/y= 44,822 
Nominal axial flow Reynolds number, Re a =Q/nDy = 48,972 
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Figure 2. (a) Description of hydrostatic recess with angled injection 

(b) Assumed pressure field within hydrostatic recess 
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eccentricity ratio e/c eccentricity ratio e/c 










eX/C 



5 recess water HJB - radial injection do:2.49mm 



eV/C 


Figure 4. Journal center locus for water - 5 recess hybrid bearing. 
Numerical predictions for increasing loads. 

(a) top - 45 degree angled injection (b) bottom - radial injection 
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centerline pressure (P-Pa)/(Ps-Pa) centerline pressure (P-Pa)/(Ps-Pa) 


5 recess water HJB - 45 deg angled injection do:2.49mm 



5 recess water HJB - radial injection do:2.49mni 



Figure 5. Numerical predictions for centerline pressure at zero load, 
Ps = 7 MPa, water - 5 recess hybrid bearing. 

(a) top - 45 degree angled Injection (b) bottom - radial Injection 


32 





centerline pressure (P-Pa)/(Ps-Pa) centerline pressure <P-Pa)/(Ps-Pa) 
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5 recess water HJB - 45 deg angled injection do:2.49mm 
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5 recess water HJB - radial injection do:2.49mm 



Figure 7. Centered recess pressure ratio (Pr-Ps) / (Ps-Pa) for water - 
5 recess hybrid bearing. Comparison to experimental results 
(a) top - 45 degree angled injection (b) bottom - radial Injection 
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WFR - whirl frequency ratio WFR * whirl frequency ratio 


5 recess water HJB - 45 deg angled injection do:2.49mm 




centered HJB : rotor speed (kcpm) 

Figure 10. Whirl frequency ratio vs. journal speed for water - 5 recess 
hybrid bearing. Comparison to experimental results. 

(a) top - 45 degree angled injection (b) bottom - radial injection 









Cross Stiffness coefficients Kxy & -Kyx [MN/m] Cross Stiffness coefficients Kxy & -Kyx [MN/m] 


5 recess water HJB - 45 deg angled injection do:2.49mm 
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Figure 11. Cross-stiffness coefficients (Kxy, -kyx) vs. journal speed 

for water - 5 recess hybrid bearing. Comparison to experimental 

results . 

(a) top - 45 degree angled Injection (b) bottom - radial injection 
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5 recess water HJB - 45 deg angled injection do:2.49mm 



5 recess water HJB - radial injection do:2.49mm 



Figure 13 . Direct damping coefficients (Cxx, Cyy) vs. journal speed 

for water - 5 recess hybrid bearing. Comparison to experimental 
results. 

(a) top - 45 degree angled Injection (b) bottom - radial injection 
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Cross damping coefficients Cxy & -Cyx [KNs/m] Cross damping coefficients Cxy & -Cyx [KNs/m] 


5 recess water HJB - 45 deg angled injection do:2.49mm 



centered HJB : rotor speed (kcpm) 


5 recess water HJB - radial injection do:2.49mm 



centered HJB : rotor speed (kcpm) 

Figure 14 . Cross-damping coefficients (Cxy, -cyx) vs. journal speed 

for water - 5 recess hybrid bearing. Comparison to experimental 
results . 

(a) top - 45 degree angled Injection (b) bottom - radial injection 





Direct inertia coefficients Mxx & Myy [kg] Direct inertia coefficients Mxx & Myy [kg] 


5 recess water HJB - 45 deg angled injection do:2.49mm 



5 recess water HJB - radial injection do:2.49mm 



Figure 15. Direct inertia coefficients (mxx, myy) vs. journal speed 

for water - 5 recess hybrid bearing. Comparison to experimental 
results . 

(a) top - 45 degree angled Injection (b) bottom - radial Injection 
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Cross inertia coefficients -Mxy & Myx[MN/m] Cross inertia coefficients -Mxy & Myx[MN/m] 


5 recess water HJB - 45 deg angled injection do:2.49mm 







WFR & Recess Pressure ratio (PR-Pa)/(Ps-Pa) 


6 recess - centered L02 HJB - angled injection do:2.33 mm 
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Figure 17. Flow rate, whirl frequency and recess pressure ratios vs. 

angle of fluid injection for 6 recess L02 hybrid bearing 
{Ps-Pa-17.9 MPa, Ts-90 K, 25 krpn) 
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6 recess - centered L02 HJB - angled injection do:2.33 mm 
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Figure 18. Shear torgue, maximum and exit temperature differences vs, 
angle of fluid injection for 6 recess L02 hybrid bearing 

(Ps-Pa-17.9 MPa, Ts-90 K, 25 krpm) 
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6 recess - centered L02 HJB - angled injection do:2.33 mm 
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Figure 19 . Stiffness coefficients (Kix=Kyy) , (Kxy— Kyx) vs. 

angle of fluid injection for 6 recess L02 hybrid bearing 

{Ps-Pa-17.9 MPa, Ts-90 K, 25 krpm) 
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Figure 20 . Damping coefficients (Cxx*Cyy) , (Cxy— Cyx) vs. 

angle of fluid injection for 6 recess L02 hybrid bearing 
{Ps-Pa-17.9 MPa, Ts«90 K, 25 krpm) 










